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We develop a trap-size scaling theory for trapped particle systems at quantum transitions. As 
a theoretical laboratory, we consider a quantum XY chain in an external transverse field acting as 
a trap for the spinless fermions of its quadratic Hamiltonian representation. We discuss trap-size 
scaling at the Mott insulator to superfiuid transition in the Bose-Hubbard model. We present exact 
f—*) ■ and accurate numerical results for the XY chain and for the low-density Mott transition in the 

' hard-core limit of the one-dimensional Bose-Hubbard model. Our results are relevant for systems 

of cold atomic eases in optical lattices. 

CO 

I. INTRODUCTION 

The achievement of Bose-Einstein condensation in diluted atomic vapors of 87 Rb and 23 Na [1] and the impressive 
progress in the experimental manipulation of cold atoms in optical lattices (see, e.g., Ref. [2] and references therein) 
have provided a great opportunity to investigate the interplay between quantum and statistical behaviors in particle 
systems. In these systems, phase transitions are phenomena of great interest, see, e.g., Refs. [3-9]. 

Phase transitions related to the formation of the Bose-Einstein condensation in interacting Bose gases at a nonzero 
temperature, as the one reported in Ref. [3], are essentially driven by thermal fluctuations, giving rise to a classical 
(_J . critical behavior, see, e.g., Ref. [10]. Quantum fluctuations play a dominant role at T = transitions, where the 
low-energy properties show a quantum critical behavior with a peculiar interplay between quantum and thermal 
fluctuations at low T, see, e.g., Ref. [11]. Quantum Mott insulator to superfiuid transitions have been observed in 
experiments with ultracold atomic gases loaded in optical lattices [4-9] . 

A common feature of the above-mentioned experimental realizations is the presence of a trapping potential V(r) 
coupled to the particle density. For example, atomic gases loaded in optical lattices are generally described by the 



T3 

q ; Bose-Hubbard (BH) model [12] 

CO '■ H BH = -5 + h h) + Z> + V ^))ni + U ni (n t - 1)], (1) 

> ' (ij) 

o . 

where (ij) is the set of nearest-neighbor sites and rij = bj&j is the particle density operator. Far from the origin the 
, potential V(r) diverges, therefore (rij) vanishes and the particles are trapped. The inhomogeneity due to the trapping 
potential strongly affects the phenomenology of quantum transitions in homogeneous systems. Correlation functions 
are not expected to develop a diverging length scale in the presence of a trap. Therefore, a theoretical description of 
how critical correlations develop in systems subjected to confining potentials is of great importance for experimental 
| investigations. 

We consider the trapping power-law potential 



k^; V(r)=vP\r\*> = (\m p , (2) 

H 

where v and p are positive constants and I = 1/v is the trap size, coupled to the particle number. Experiments are 
usually set up with a harmonic potential, i.e., p = 2. 

Let us consider the case in which the system parameters are tuned to values corresponding to the critical regime of 
the unconfined system. In the presence of a confining potential, the critical behavior of the unconfined homogeneous 
system can be observed around the middle of the trap only in a window where the length scale £ of the critical modes 
is much smaller than the trap size, but sufficiently large to show the universal scaling behavior. If £ is large, but not 
much smaller than the trap size, the critical behavior gets somehow distorted by the trap, although it may still show 
universal effects controlled by the universality class of the phase transition of the unconfined system. In this paper 
we investigate this issue within quantum transitions. 

In Ref. [13] we have shown that the critical behavior of trapped systems at classical continuous transitions can be 
cast in the form of a trap-size scaling (TSS), resembling the finite-size scaling theory for homogeneous systems [14, 15], 
but characterized by a further nontrivial trap critical exponent 9, which describes how the length scale £ depends on 



2 



the trap size at criticality, i.e., £ ~ l e , and which can be estimated using renormalization-group (RG) arguments. TSS 
was supported by numerical results for some lattice gas models. 

In the present paper, we extend the study of the effects of trapping potentials to quantum critical behaviors at T = 
transitions. We put forward a TSS scenario to describe how critical correlations develop in large traps, analogous 
to the one outlined in Ref. [13] for classical continuous transitions. We then check the validity of the TSS scenario 
within quantum transitions by determining the trap-size dependence in a few quantum models. As a theoretical 
laboratory, we consider a quantum XY chain in an external space-dependent transverse field, acting as a trap for 
the spinless fermions of its quadratic Hamiltonian representation. We present exact and very accurate numerical 
results which fully support TSS. Moreover, we discuss TSS at the Mott insulator to superfluid transition in the BH 
model. In particular, we present analytical results for the low-density Mott transition in the hard-core limit of the 
one-dimensional BH model. The results again support TSS, although the corresponding TSS functions show peculiar 
behaviors, like discontinuities in the scaling particle density, which are related to the quantum nature of the transition. 

The paper is organized as follows. In Sec. II we outline the main ideas of the TSS theory. In Sec. Ill we study the 
XY chain in a space-dependent transverse field, which gives rise to an inhomogeneity analogous to the one arising 
from a trapping potential in particle systems. We show how TSS emerges by analytical and very accurate numerical 
results. In Sec. IV we consider the BH model (1) in a confining potential. We discuss the general features of the 
Mott insulator to superfluid transitions, and present the RG arguments to derive the corresponding trap exponent 9. 
We then present results for the TSS at the low-density Mott transition in the one-dimensional hard-core BH model. 
Finally, in Sec. V we draw our conclusions. 



II. TRAP-SIZE SCALING AT QUANTUM TRANSITIONS 

Quantum phase transitions arise from a nonanaliticity of the ground-state energy, where the gap A = E\ — Eq 
vanishes. In addition, continuous quantum transitions have a diverging length scale £, which gives rise to peculiar 
scaling properties, see, e.g., Ref. [11]. 

In a standard general scenario, a quantum T = transition of a homogeneous d-dimensional system is generally 
characterized by one relevant parameter /x, with critical value fi c and RG dimension = 1/v, and a given dynamic 
critical exponent z. Interesting physical examples are the Mott insulator to superfluid transitions in BH models. The 
quantum critical behavior is characterized by a vanishing energy scale, for example the gap A, and a diverging length 
scale £, behaving respectively as A <~ \p,\ zu and £ <~ Ifll^", where p, — \i — fj, c , thus A ~ £~ z . 

The confining potential gives rise to a space inhomogeneity, thus changing the scaling behavior of the homogeneous 
system, which could be only recovered in the limit of large trap size (keeping fixed the other parameters of the system). 
We want to describe how the critical correlations develop in a trap. Our starting point is a scaling Ansatz which 
extends the scaling laws of homogeneous systems at quantum transitions, to allow for the presence of a confining 
potential like (2). We write the scaling law of the singular part of the free energy density at the quantum transition 
as 

F(n,T,v,x) = F(fib y »,Tb z ,vb y \x/b), (3) 

where b is any positive number, y v is the RG dimension of the trap parameter v, cf. Eq. (2), and x is the distance 
from the middle of the trap. We are neglecting irrelevant scaling fields, because they do not affect the asymptotic 
behaviors. Then, fixing vb Vv = 1 and introducing the trap size I = v , we obtain the following trap-size scaling 
(TSS) 

F = r e V+*)T(pl e l v , Tl 0z , xl~ e ), (4) 

where v = 1 jy^ and 9 = l/y v is the trap exponent. 9 depends on the universality class of the quantum transition. It 
can be computed by evaluating the RG dimensions of the corresponding perturbation. 

Analogously, one can derive the TSS of other observables. Any low-energy scale at T = 0, and in particular the 
gap, is expected to behave as 

A = r e *V(jil e ' v ), (5) 

with T>(y) ~ y zv for y — > to match the scaling A ~ \± zv in the absence of the trap. The correlation length £ around 
the middle of the trap, or any generic length scale associated with the critical modes, behaves as 

t = l e X{jll 9/v ,Tl 9z ), (6) 
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where X(y, 0) ~ y~ v for y — > 0. This implies that at the T — quantum critical point the trap induces a finite length 
scale: £ <~ l e . In the case of a generic local operator O(x), with RG dimension y , we expect that its expectation 
value and equal-time correlator behave as 

(0(x)) = r Bv °o{pi e/ \ ti 6z , x r e ), (7) 
(O(x)0(o)) c = r 2e y°g {pi e,v , Tf\ xl- e ), (8) 

where x measures the distance from the origin, i.e., the middle of the trap. In the above scaling formulae we have 
neglected scaling corrections due to irrelevant perturbations, and possible analytic contributions. 

When p — > oo the effect of the trapping potential is equivalent to confine a homogeneous system in a box of size 
L = 21 with open boundary conditions, whose behavior can be described by a standard finite-size scaling [10, 14]. 
Therefore we must have that 9 — > 1 when p — > oo. 

The TSS theory outlined in this section provides an effective theoretical framework to describe quantum critical 
behaviors in confined systems. However, it is important to check the validity of its scaling Ansatz, because the quantum 
nature of the phenomenon may lead to subtle effects. For this purpose, in the following sections, we study the trap-size 
dependence at some specific quantum transitions of one-dimensional systems in the presence of a space-dependent 
confining potential. 



III. TRAP-SIZE SCALING IN THE QUANTUM XY CHAIN 

The quantum XY chain in a transverse field is a standard theoretical laboratory for issues related to quantum 
transitions, see, e.g., Ref. [11]. In this model, an inhomogeneity analogous to the one arising from a trapping potential 
in particle systems can be achieved by considering a space-dependent transverse field, i.e., 

H ™ = E + -y)°i°i+i + (! - 7)*M+i] - 1*1 V( Xi )at, (9) 

i 

where o % are the Pauli matrices, < 7 < 1 and V(x) — v p \x\ p = (\x\/l) p , where v and p are positive constants 
and I = v^ 1 is its length scale. This model can be mapped into a quadratic Hamiltonian of spinless fermions by a 
Jordan- Wigner transformation, 

H = ^[ctAiM + \{c\B ijC ) + H.c.)], 

Aij = 2Sij — Si+l,j — + 2Q(%i)o~ij, 

B ij = -7 - kj+i) > ( 10 ) 

Q(x) = ft + V(x), p, = /j,-l. 

In this picture fx plays the role of chemical potential for the c-particles, and the space-dependent field V(x) acts as a 
trap for the c-particles, making their local density (rij) = (qcj) vanish at large distance. In the following, I = w _1 will 
be considered the trap size. 

In the absence of the trap, the model undergoes a quantum transition at p, = in the two-dimensional Ising 
universality class, separating a quantum paramagnetic phase for p, > from a quantum ferromagnetic phase for 
p, < 0. Around p — 0, the quantum critical behavior shows a diverging length scale, £ ~ \p\~ u , and a vanishing energy 
scale, A ~ , where z and v are universal critical exponents: z = 1 and v = l/y M = 1 (y^ is the RG dimension of 

m)- 

In the presence of the confining potential, the critical behavior can be observed only in the limit of large trap size. 
The trap exponent 8 controlling such limit can be determined by analyzing the RG properties of the corresponding 
perturbation at the critical point, which can be represented by Py = j d d xV(x)(f>(x) 2 , where <j>(x) is the order- 
parameter field of the <\> theory which describes the behavior of the critical modes [10, 13]. Since the RG dimensions 
of the potential V(x) and the energy operator <f) 2 , respectively yy = py v — p and y^i = d + z — y^, are related by 
yv + = d + z, we obtain py v — p = y^, and therefore, since z = 1 and y^ = 1, 

6=l/y v =p/(p+l). (11) 

Notice that 9 — > 1 when p — > 00; this was to be expected, since when p — > 00 the system becomes equivalent to a 
homogeneous chain with —l<x<l and open boundary conditions. 

We mention that some issues related to the presence of a space-dependent transverse field in the quantum XY chain 
model have been addressed in Refs. [16-19]. 
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A. The trap-size scaling limit 



The Hamiltonian (10) can be solved by exact numerical diagonalization, even in the presence of the trapping 
potential, following, for example, Ref. [20]. We look for new canonical fermionic variables r\ k = 9kic\ + hkiCi which 
diagonalize the Hamiltonian, i.e., such that 

H = UkVkVk, u k > 0, (12) 
k 

and the ground state satisfies r?fc|0) = for any k. The diagonalization is achieved by introducing two sets of 
orthonormal vectors <p k and ij) k with components <j) k i = Ski + h k i and Vfei = 9ki — h k % respectively, satisfying the 
equations 

(A + B)4> k = cjfcVfc, (A- B)ip k =uj k (j) k , (13) 
where the matrices A and B have been defined in Eq. (10). Thus, we have 

(A-B)(A + B)<p k =uU k . (14) 
Once the vectors 4> k have been computed, one can easily determine tp k and reconstruct the original c-operators as 

){4>ki + 4>ki)Vk + {<t>ki - ipki)%, (15) 
k 

from which one can compute the correlation functions of the c-operators, by straightforward calculations. 

In order to show how TSS emerges, we consider the continuum limit of Eq. (14). A similar approach was developed 
in Ref. [16] to study the critical behavior of the quantum XY chain in the presence of a gradient perturbation. By 
rewriting the discrete differences in terms of derivatives, near the critical point p, — and for sufficiently small values 
of k (this is required by the smoothness hypothesis underlying the continuum limit), we obtain 

[4Q(xf + 4^d x Q(x) - 4j 2 2 x - 4Q{x)d 2 x - 2d x Q(x)d x - 2d 2 x Q(x) + ...] <p k {x) = w 2 k <i> k {x). (16) 

This equation has a nontrivial TSS limit: by rescaling 

x = 7 V(i+p) F /(i+p) X , 

/ u = 7 p/( 1 +p)/-p/( 1 +f) Mr , (17) 

cj k = 2 7 p/( 1+ p)r p /( 1+p ^ fe , 

and keeping only the leading terms in the large-l limit (and for small u) k ), we obtain 

( Mr + x p - d x ) ^ + x p + d x ) MX) - nlM*)- ( 18 ) 

Analogously, we obtain the equation for the function ip k (x), i.e., 

{Hr + X"-dx)MX) = n k </> k (X). (19) 

It is worth noting that the dependence on 7 in Eqs. (18) and (19) disappears, showing the universality of TSS with 
respect to changes of the values of 7 (of course, excluding the value 7 = 0). One can easily check that the next- to- 
leading terms in the large-trap limit give rise to 0(l~ e ) scaling corrections. Corrections due to irrelevant perturbations 
already present in the homogeneous system are expected to be more suppressed: the exponent associated with the 
leading scaling corrections in generic systems belonging to the two-dimensional Ising universality class is w = 2 [10, 21], 
thus leading to 0{l~ 2e ) corrections in TSS. 

The solution of Eqs. (18) and (19) allows us to determine the critical correlations of the original c-operators, 
using Eq. (15). One can easily determine the large-distance behavior of the functions (j) k {x) and ip k (x): they decay 
exponentially as <~ exp(— a|x| p+1 ) apart from prefactors. Their large-distance decay determines the large-distance 
behavior of the scaling correlation functions of the c operators. 

A natural question concerns the robustness of the TSS results with respect to variation of the trapping potential 
from a simple power law. For example, let us consider a potential given by 

u{x) = {\x\/iy + b{\x\/iy, q> P , (20) 

where b > is a constant. A simple analysis based on the rescalings (17) shows that the TSS limit is determined by 
the smallest power law, while the higher power gives rise to 0(/ _ ( < ? _ p)/( 1 +p)) scaling corrections. 
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B. Trap-size scaling of observables 



In the following we report results obtained by exact diagonalization at fixed trap size, and compare them with the 
TSS derived in Sec. II. The exact numerical diagonalization in the presence of the trap is done on a chain of length 
L with open boundary conditions, with L large enough to have negligible finite-size effects. 



1. Trap-size scaling of low-energy scales 




FIG. 1: (Color online) TSS of the energy differences Ai = E\ — Eq and A2 = E2 — Eo for p = 2 (above) and p — 4 (below) vs. 
/i r = 7~ s i e /x, for several values of 7 and for I > 10. Numerical diagonalization results clearly approach universal TSS functions 
in the large-Z limit (represented by full lines and obtained by extrapolations), with 0{l~°) scaling corrections (larger at small 
7 and for higher levels). 



The universal TSS limit obtained after the rescalings (17) implies the following asymptotic behavior for any low- 



6 



energy scale: 

A« 7 e r e 2?( Mr ), Mr = 7 -^A, (21) 

which is approached with 0(l~ e ) scaling corrections. This proves the scaling behavior (5) obtained by RG arguments. 

In Fig. 1 we plot the differences among the first few energy levels, Ai = E\ — Eo and A2 = E2 — Eo, obtained 
by numerical diagonalization, for trapping potentials with p = 2 and p = 4 and for several values of [i and 7. The 
numerical results clearly approaches a universal TSS behavior in the large-/ limit, as predicted by Eq. (21). The 
large- 1 [i r | behavior of the scaling functions T>i related to A, are T>i <~ 2u r for [i r — > +00, V\ — > for /j, r — > —00, 
and T>2 ~ 2\/j, r \ for /U r — > ±00; they are consistent with the behaviors of the corresponding quantities in the absence 
of the trap. The vanishing of T>\ for fj, r — > —00 is related to the degeneracy of the ground state in the quantum 
ferromagnetic phase. 



2. Trap-size scaling of the particle density and its correlators 

Other interesting quantities are the particle density and its correlation function. At the quantum critical point, 
the particle density is expected to behave analogously to the energy density in the standard Ising model [13], which 
presents leading analytic contributions arising from the analytic part of the free energy. At p, = and in the middle 
of the trap, we expect 

(n ) = (c c ) = p c ( 7 ) + a nl e /n~ 6 + ..., (22) 

where we used the fact that the RG dimension of the particle density operator is y n = d+z — y^ = 1, and the rescalings 
(17) to guess the dependence on 7 of the amplitude of the scaling term. a n is a constant that does not depend on 
7. p c {l) is the critical value of the homogeneous system without trap, for example p c (l) = 0.18169011... [22]. Our 
numerical results for several values of 7 support the scaling behavior (22). Scaling corrections are (consistent with) 
integer powers of l~ e . 

We have also studied the dependence on the distance from the middle of the trap, which must be an even analytic 
function of x at fixed trap size, but it may present a nontrivial scaling for large I. Analogously to the scaling at the 
origin, we expect it to be given by a sum of an analytic behavior and a nonanalytic term in the trap size, such as 

(n x - n ) = p ns (x/l) + a n r K(xr ), (23) 

where p ns is a nonsingular term which we expect to be a function of x/l. Let us compare this Ansatz with the results 
from exact diagonalization. At small x, the latter are consistent with the expansion 

(n x - n ) « b 2 (l)x 2 + 6 4 (0z 4 + (24) 

with 

6 2 (0 = b 21 r 2 ini + b 22 r 2 + ..., 6 4 (0 = &4ir 10/3 + ... (25) 

for p = 2 (dots indicate terms that are more suppressed in the large-Z limit by power laws), and 

62(0 = b 21 r 12/5 + b 22 r 16/5 + ..., 6 4 (0 = &4ir 4 in/ + 6 42 r 4 + ... (26) 

for p = 4. These results should be compared with the expansion of Eq. (23) in powers of x. They are in substantial 
agreement; the only difference is given by the presence of logarithms. However, this should not be surprising, because 
logarithmic corrections to power laws are peculiar of the two-dimensional Ising universality class, see, e.g., the behavior 
of the specific heat, and arise [23] from integer resonances of the RG dimensions of different perturbations, see also 
Ref. [24]. Note that the logarithms appearing in the expansion (24) show up in the x 2 term for p = 2 and in the 
x 4 term for p = 4, where the trap-size power-law expansions of the two terms in Eq. (23) present an accidental 
degeneracy, just because 39 = 2 for p = 2 and 59 = 4 for p = 4. The large-distance behavior of the particle density 
turns out to be dominated by the nonuniversal analytic term, indeed at large distance we find 

K>*4(f)" 2P > (27) 

while the scaling nonanalytic part is expected to be exponentially suppressed, as suggested by the large-distance 
behavior of the solutions of Eq. (18). 
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FIG. 2: (Color online) TSS of G„{x), cf. Eq. (28), at p, = vs. X = ^ 0/p r e x, for p = 2 (above) and p = 4 (below). Numerical 
diagonalization results clearly approach universal TSS functions in the large-Z limit (represented by full lines and obtained by 
extrapolations) . 



The static particle-density correlator is not affected by analytic backgrounds; therefore, according to Eq. (8), we 
expect that at /2 = and for x ^ 

G n (x) = (n Q n x ) - (n )(n x ) » j 2e/p r 2e g n (X), (28) 

where Q n (X) is a universal function. This is confirmed by the results of numerical diagonalization, as shown in Fig. 2 
for p — 2 and p = 4. At small X, Q n (X) ~ l/^ 2 , which is the behavior in the absence of trap, while at large X it 
decays very rapidly. 

In the limit p — ¥ oo the system becomes equivalent to a homogeneous finite-size chain with —l<x<l and open 
boundary conditions. Thus, in this limit we can exploit conformal field theory (CFT), by mapping the half-plane 
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result for the two-dimensional Ising universality class [25] into the strip, obtaining the correlation function [26] 

_ , . 7r 2 cos(7rX/2) „ ,, . „. 

G ^ x) = ^^My X = X,L (29) 

Since 9 —> 1 for p —> oo, this result agrees with TSS. Of course, the corresponding scaling function Q n (X) = l 2 G n (x) 
is defined for \X\ < 1, and vanishes for X — ±1. It is shown in Fig. 3. 



3. Trap-size scaling of two-point correlation function 



We have also computed the static correlation function of the spin operator a x at T — and /j, c , following essentially 
the method of Ref. [27]. According to Eq. (8), G s (x) behaves as 

G s (x) = (o%<4) = a s l- e ^g s (X), (30) 

where we used the fact that the RG dimension of the spin operator a x is rj/2 with rj = 1/4. The constant a s is not 
universal, it is expected to depend on 7. Our results from numerical diagonalization are perfectly consistent with 
Eq. (30): we find that, introducing a 7- (and p-) dependent normalization c, determined phenomenologically, the 
combination cl 8ri G s (x) plotted vs. X = 7~ 1 /( 1 +p)/~ 6 'x, for different values of x (not too small, e.g., x > 4), I (large), 
and 7, falls on a single curve; see Fig. 4. Q S (X) diverges as X^ 1 / 4 for X — > (corresponding to the behavior in the 
absence of the trap) and decay rapidly at large X, as exp(— ax^ p+1 ^), in agreement with the large distance behavior 
of the solutions of Eq. (18). 

One may also consider the susceptibility and second moment correlation length defined as 



X 



dr G s (r) » c x l {d ~i )e + B = a x P/ 4 + B, (31) 



e = y \ drr 2 G s (r)^a e l 29 , (32) 

where B is a background constant term already present in the homogeneous system without trap [10]. Since the 
normalization of G s (x) cancels out, £ must scale as x, i.e., 

Z = a a e /n e [l + 0(r% (33) 



with ci£ dependent on p but not on 7. The above scaling behaviors have been verified by numerical computations. 
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FIG. 4: (Color online) TSS of the two-point function G s (x), cf. Eq. (30), for p = 2 (above) and p — 4 (below). The universal 
TSS functions are represented by full lines and obtained by extrapolations in the large-Z limit. 



Finally, the p = oo limit of the correlation function (30) can be computed using CFT, analogously to the density- 
density correlation function (29), obtaining [26] 



G s (x) 



-,1/8 



Z 2 16F 2 sm 2 (nX/4) 



y/T+Y- VT^Y 



s/2 



Y = Jl -tan 2 (7rX/4), X = x/l. 
The corresponding scaling curve g s (X) = l^ 4 G s {x) is plotted in Fig. 5. 



(34) 



4- Trap-size scaling of quantum entanglement 
We now discuss quantum entanglement in the presence of the confining potential. 
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FIG. 5: (Color online) The scaling function Q S {X) = l 0r >G s (x) vs. X = x/l B in the limit p ->■ oo where 9 -> 1, cf. Eq. (34). 



We divide the chain in two parts of length and L — 1^. At /i = 0, in the absence of the trap and for open 
boundary conditions, the entanglement entropy is given by [28] 

S(l A ;L) ^ -In [Lsm(nl A / L)}+ E (35) 
6 

where c = 1/2 is the central charge. The constant E depends on 7: 

E= l + iV n(27A) (36) 

where a = 0.478558... [29]. 

In the presence of a trapping potential, and for sufficiently large lattice sizes L, we obtain 

5(L/2;L)« ~ln£ e + £, (37) 

D 

which defines an entanglement length £ e at the critical point. Note that E is the same constant appearing in Eq. (36). 
The ratio between £ e and the correlation length £ defined from the correlation function of a x , cf. Eqs. (32) and (33), 
should be universal in the large-Z limit, thus independent of 7. This implies that at criticality £ e has the asymptotic 
behavior 

£ e = a e7 e /n e [l + 0(r% (38) 

with a e dependent on p but not on 7. 

The above scenario is fully supported by numerical computations, using the technique of Refs. [18, 30]. 

C. Trap-size dependence for /1 < 1 

Finally, it is worth discussing the trap-size dependence at values of /1 < 1, i.e., /2 < 0, within the quantum ordered 
phase. In this case the TSS theory outlined in Sees. II and III A does not apply, because the system is not critical 
in the middle of the trap. Around the middle of the trap we have a quantum ferromagnetic phase, up to the point 
ft + V{x) = 0, i.e., up to \x c \ — {~{jL) 1 / p l, around which critical fluctuations develop. Then, as before, at larger 
distances the system is in the quantum paramagnetic phase and the particle density vanishes for |sc| — > 00. 

In Fig. 6 we plot the differences of the lowest energy levels, i.e., Aj = E\ — E and A 2 = E 2 — E , for a few values of 
fi < 1 for the harmonic potential, i.e., p — 2. Ai vanishes rapidly, apparently as Ai ~ exp(— nl), where the constant 
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k depends on /j,. This behavior is clearly related to the degeneracy of the ground state in the quantum ferromagnetic 
phase in the absence of trapping potential. 

It is interesting to note that also A2 appears to vanish in the large-? limit, although at a much slower rate. Indeed, 
our data for fi not too close to the critical value fi = 1 are consistent with the power-law behavior A 2 <~ Z -1 / 2 . The 
vanishing of A 2 should be related to the developing of critical fluctuations at the points x c = ±(—fi) 1 / p l where the 
sum p, + V(x) vanishes. Since around x c the potential V(x) = \x\ p /l p can be approximate by a linear potential, i.e., 

V(x) nu(x-x c ), u = p(-fi) 1/p /l, (39) 

the exponent 1/2 should be related to the RG scaling dimensions of the parameter u of the linear potential. Using 
RG arguments analogous to those leading to Eq. (11), we obtain y u — 2, see also Ref. [16]. This explanation is further 
supported by the fact that the data for A2 for different values of n < 1 appear to merge to a unique curve as a 
function of u for sufficiently large trap size I, as shown in Fig. 6 (fi — 0.99 clearly displays crossover effects because it 
is very close to the critical value \l = 1). 



IV. TRAP-SIZE SCALING IN THE BOSE-HUBBARD MODEL 



A. General features 



We now discuss TSS within the Bose-Hubbard (BH) model (1) at the Mott insulator to superfluid transitions. 
A review of results for the BH model can be found in Ref. [11]. In the homogeneous BH model without trap, 
the low-energy properties of the transitions driven by the chemical potential \i are described by a nonrelativistic 
U(l)-symmetric bosonic field theory [31], whose partition function is given by 



= y [L>0] exp ^- ^ 7 dtd d xc}j, 



(40) 



C c = cj>*d t cj> + ^|V0| 2 + r\4>\ 2 + u\^\\ (41) 

where r ~ /j — /j, c . The upper critical dimension of this bosonic theory is d = 2. Thus its critical behavior is mean field 
for d > 2. For d = 2 the field theory is essentially free (apart from logarithmic corrections), thus the dynamic critical 
exponent is z = 2 and the RG dimension of the coupling \i is — 2. In d = 1 the theory turns out to be equivalent 
to a free field theory of nonrelativistic spinless fermions, from which one infers the RG exponents z = 2 and y^ = 2. 

The special transitions at fixed integer density (i.e., fixed ji) belong to a different universality class, described by a 
relativistic U(l)-symmetric bosonic field theory [31], whose Lagrangian is 

£ c = |9 t 0| 2 + « 2 |V0| 2 +r|^| 2 + U |0| 4 , (42) 

for which z = 1 and y^ — 1/vxy where vxy is the correlation length exponent of the D = d+ 1 XY universality class. 
Thus vxy = 1/2 for d = 3 (i.e., mean-field behavior apart from logarithms), vxy = 0.6717(1) in the d = 2 case [32], 
and formally v = 00 for the Kostcrlitz-Thoulcss transition [33] at d = 1. 

In the presence of a confining potential, theoretical and experimental results have shown the coexistence of Mott 
insulator and superfluid regions when varying the total occupancy of the lattice, see, e.g., Refs. [8, 12, 34-37]. However, 
at fixed trap size, the system does not develop a critical behavior with diverging length scale [34, 35]. Criticality is 
recovered only in the limit of large trap size, where the critical behavior can be described in the framework of the 
TSS theory. The trap exponent 9 can be determined by an analysis of the corresponding RG perturbation, i.e., 



d d xdtV{x)\(f)(x)\ 2 , (43) 
leading to the relation 

PVv -P = d + z - y W 2 = Vfl , (44) 
where y^ is the RG dimension of the coupling r of the quadratic term of Lagrangian (41) or (42). Therefore, 

* S 1 = _E_. (45) 

Vv p + y^ 



By replacing the corresponding value of y^ , this relation yields the value of 9 for each specific transition. 
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FIG. 6: (Color online) Trap-size dependence at /j, < 1, i.e., in the quantum ferromagnetic phase, of the difference of the lowest 
energy levels for n < 1 and p = 2: A x = - E vs. i (above) and A 2 = E 2 - E vs. u 1/2 = 2 1/2 (-p,) 1/4 r 1/2 (below). 



B. TSS at the low-density Mott transition in the ID hard-core Bose-Hubbard model 

We now show how TSS emerges at the low-density Mott transition of the one-dimensional BH model, which is also 
of experimental relevance in optical lattices, see, e.g., Refs. [2, 5, 6, 9]. In particular, we consider its hard-core limit 
U — > oo, which implies that the particle number is restricted to the values n% — 0, 1. In this limit, analytical results 
can be obtained by exploiting exact mappings into the so-called XX chain model in the presence of a space-dependent 
external field, 

i i i 
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(which is the Hamiltonian (9) with 7 = 0). The Pauli spin matrices are related to the boson operators of the Bosc- 
Hubbard Hamiltonian i?BH, cf. Eq. (1), by af = b\ + bi, of = i(b\ — of = 1 — 2&jtV Actually we have that 
Hxx = %Hbh(J = l,f — > 00). One can then map it into a model of free spinless fermions by a Jordan- Wigner 
transformation, given by Eq. (10) for 7 = 0, see, e.g., Ref. [11]. 

In the absence of the trap, the ID hard-core BH model has three phases: two Mott insulator phases, for /x > 1 with 
(rii) = and for fj, < — 1 with (rij) = 1, separated by a gapless superfluid phase for < 1. Therefore, there are two 
quantum transitions at \i = ±1, with z = 2 and = 1/v = 2. In the following we consider the low-density Mott 
transition at it = 1, for which we can analytically show the existence of the TSS limit, and compute the corresponding 
scaling functions. The results provide a further analytical check of the TSS theory in quantum transitions, in a case 
with dynamic exponent z ^ 1, whose quantum critical behavior is not described by a d + 1-dimensional conformal 
field theory. 

The trap-size dependence, and corresponding TSS, turns out to be more subtle at the (n») = 1 Mott insulator to 
superfluid transition, i.e., at \i = — 1. This is essentially due to the fact that at fj, = — 1 there is an infinite number of 
level crossings as I — » 00. 1 Results will be presented elsewhere. 

In the fermion representation the Hamiltonian can be easily diagonalized: introducing new canonical fcrmionic 
variables r\ k = tfikiCi, k = 0, 1, . . . , where ip k i satisfies 

Aij(p kj = ojkfkj (47) 

and Aij is the matrix defined in Eq. (10), we obtain 

H = J2^kvlvk- (48) 

k 

The ground state contains all 77-fermions with ojk < 0, therefore the gap is 

A = min|w fe |. (49) 

k 

A nontrivial TSS limit around p, = fj, — 1 = 0, i.e., at the transition between a low-density superfluid and the empty 
vacuum state (named (rij) = Mott phase above), is obtained (for small \tJk\) by rescaling 

x = ip/<*+p) x , ^ = r 2 ^ 2 +*V, u; k = r 2 p/^n k . (50) 

Indeed, neglecting terms which are suppressed in the large-? limit, Eq. (47) becomes 

(2XP - 2 x ) Vk {X) = (fi fc - 2^r)MX)- (51) 

This shows that 6 = p/(2 + p), in agreement with Eq. (45). 

Moreover, this implies that any energy scale, and in particular the gap A = E\ — E , must behave as 

A w l~ 2e V(p, r ), fir^l 28 ^ (52) 

which agrees with the RG scaling equation (5), since z — 2 and is = 1/2. For p = 2, by solving Eq. (51), we obtain 

V(fi r ) = min|2 3/2 (/c + l/2) + 2ii r |, fc = 0, 1, . . . ; (53) 

k 

T>(/j, r ) is a triangle wave for /i r < and it is linear for /j, r > —l/y/2. 

The TSS of the particle density in the middle of the trap is obtained by computing 

<no>=r'5>*(0)| 2 <^fc>, (54) 
k 



For fi < 1, the ground state contains all the r)-fermions diagonalizing the Hamiltonian with u) k < 0, cf. Eq. (48). In the presence of the 
trapping potential (2), level crossings of the lowest states occur in the ji-l plane separating the region with N = = n from 

N = n+ 1. Since, in the absence of the trap potential and for \x < 1, the ground state has a finite density N/L > 0, for fixed pi < 1 the 
lowest states show an infinite number of level crossings as I — > oo (after L — > oo) where the gap vanishes. 
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FIG. 7: (Color online) TSS of the gap (below) and the particle density in the middle of the trap (above) for p — 2 and several 
values of I > 10 and fl. The lines correspond to Eqs. (53) and (55). Scaling corrections turn out to be very small. 



where tpk(X) are the normalized eigenfunctions of Eq. (51); = 1 if fife < and otherwise; since (fk(X) = 

(— 1) <^/-(— X), only even fcs contribute. For p = 2 we obtain the sum 

i 9 (n ) = (2 i / 4 /v^)E toy - m*ttw (55) 

over integer j > satisfying 2 1 / 2 (2j + l/2) + u r < 0. Again, this result agrees with the TSS theory, taking into account 
that the RG dimension of the particle density is y n = d+ z — = 1 . Results from numerical diagonalization at p = 2 
are shown in Fig. 7; they fully support the above TSS behaviors. Note the peculiar plateaux and the discontinuities 
in the particle density at negative values of the scaling variable p, r = l 20 Ji. For p r — > — oo, (no) ~ -\/|2/2|/7r, which 
matches the critical behavior for p, < in the absence of the trap [11]. 

Numerical results for higher power laws of the potential turn out to be in full agreement with the predictions of 
TSS and are qualitatively similar to the results for p — 2, including those for p — > oo which effectively corresponds to 
the finite-size scaling of the homogeneous system without confining potential. 



V. CONCLUSIONS 



We have developed a trap-size scaling (TSS) theory for confined particle systems at quantum transitions. We 
propose an extension of the scaling laws of the homogeneous system to allow for the confining potential in the large 
trap-size limit. The quantum critical behavior of the confined system is cast in the form of a TSS, resembling finite-size 
scaling theory [10, 14, 15], with a nontrivial trap critical exponent 9, which describes how the length scale £ of the 
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critical modes diverges with increasing trap size, i.e., £ <~ I , at the quantum critical point. This provides a theoretical 
framework to describe the quantum critical behavior of confined systems. 

We have shown by explicit analytical computation how TSS emerges in the quantum XY chain, where an inhomo- 
geneity analogous to the one arising from a trapping potential in particle systems can be achieved by considering a 
space-dependent transverse field. 

Moreover, we have discussed this issue within the Bose-Hubbard model with a confining potential, which is relevant 
for the description of cold atomic gases in optical lattices. In particular, we have presented some analytical results for 
the low-density Mott transition in the hard-core limit of the one-dimensional BH model. The results again support 
TSS, although the corresponding TSS functions show a peculiar behavior, like discontinuities in the scaling particle 
density, which are clearly related to the quantum nature of the transition. 

Helpful discussions with P. Calabrese are gratefully acknowledged. 
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